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SUMMARY  OF  PROGRESS  FOR  THE  THIRD  AND  FOURTH  QUARTERS,  1983 


This  is  a  summary  of  highlights  of  progress  made  by  our  group,  for  use  by 
DOE  and  ONR  contractors.  The  details  appear  in  our  Quarterly  Progress 
v  Report. 

>  /A  ■  SJ~0  (»  ►-  '  ,  f  4  «-  f  ■'  *  6  t  ’  *  h  •  *  /, 

Our-grour-uses- theory  and  simulation  as  tools  in  order  to  increase 
the  understanding  of  instabilities,  heating,  transport  and  other  phenomena 
in  plasmas.  ,^We  also  work  on  the  improvement  of  simulation  both  theoretically 
and  practically. 

Our  group  consists  of  Prof.  Birdsall,  Dr.  Crystal,  Dr.  S.  Kuhn  (visiting 
from  Innsbruck)  and  six  graduate  students.  We  are  aided  by  professionals 
from  the  Lawrence  Livermore  National  Lab  and  the  Space  Sciences  Lab. 

I.  Theory  and  Simulation 

A.  Simulation  of  the  Ion-Beam-Driven  Drift  Instability  in  a  Magnetic  Trap 

Complete  report  done  and  sent. 

B.  Application  of  Nonlinear  Constants  of  Motion  in  a  Single  Electromagnetic 
Wave  to  the  Study  of  the  Alfven-Ion-Cyclotron  Instability 

Report  now  complete  and  sent. 

C.  Theory  of  Nonmonotonic  Double  Layers 

Report  now  complete  and  sent. 

D.  Magnetized  Sheath  Simulations  (joint  with  IPP,  Nagoya) 

For  the  magnetic  field  normal  to  a  wall  (same  as  the  magnetic  field 
in  the  ld3v  model,  the  transient  rarefaction  wave  is  observed,  moving 
into  the  plasma  at  c  .  For  the  BQ  field  of  angle  ij;  f  0  to  the  normal, 
the  ion  drift  velocity  at  the  wan  drops  off  roughly  as  cosiJj  for  suffi¬ 
ciently  large  (uci/wpi);  the  wall  potential  is  rather  insensitive  to 

measured  after  the  initial  transient.  0<^£80°  was  used.  (In  previous 
QPR's,  we  showed  for  =  90°,  that  the  plasma  becomes  negative  with 
respect  to  the  wall). 

E.  Saturation  Characteristics  of  Counterstreaming  Warm  Electrons 

The  distribution  function  at  v  =  0,  f(0),  is  found  to  overshoot  the 
value  for  a  Maxwellian  for  vthermal/vdrift  '  °'3  an,i  lJrift  tow,rd  the 
Penrose  stable  value  for  .  >  0.5. 


F.  Electron  Bernstein  Wave  Investigations  Regarding  Linear  and  Non-Linear 
Damping 

For  off  perpendicular  propagation,  spatial  (Landau)  damping  has 
been  observed,  as  predicted. 

6.  Simulation  of  the  Classical  Pierce  Diode 

This  well-known  instability  has  been  examined  theoretically,  through 
several  eigenmodes,  which  has  been  verified  by  simulations. 

H.  Simulation  of  the  Plasma-Sheath  and  Presheath  Regions 

Comparisons  of  simulations  with  time-independent  theory  for  m./m  = 

40  and  100  for  4>wa]-|  are  excellent.  This  work  continues.  1  e 


II.  Code  Development 

A.  Plasma  Diode:  Id  Vlasov  simulation  (GASBAG  Code) 

Report  completed  and  sent. 

B.  PDW1,  Plasma  Device  Code 

Report  completed  and  sent. 

Publ ications 

Two  journal  articles  and  five  reports  were  published. 

Distributed  with  this  Progress  Report  are: 

"Linear  Longitudinal  Oscillations  in  Collisionless  Plasma  Diodes  with  Thin 
Sheaths.  Part  I.  Method,"  by  S.  Kuhn,  Physics  of  Fluids,  27^(7),  July  1984, 
pp.  1821-1833. 

"Linear  Longitudinal  Oscillations  in  Collisionless  Plasma  Diodes  with  Thin 
Sheaths.  Part  II.  Application  to  an  Extended  Pierce-Type  Problem,"  by 
S.  Kuhn,  Physics  of  Fluids,  22(7),  July  1984,  pp.  1834-1851. 

"PDWl  User's  Manual,"  by  W.  S.  Lawson,  ERL  Technical  Memorandum  No.  UCB/ERL 
M84/37,  27  April  1984. 

"Simulation  of  the  Ion-Beam-Driven  Drift  Instability  in  a  Magnetic  Trap  I," 
by  V.  A.  Thomas,  W.  M.  Nevins,  and  Y-J.  Chen,  ERL  Technical  Memorandum  No. 
UCB/ERL  M84/45,  13  June  1984. 

"Simulation  of  the  Ion-Beam-Driven  Drift  Instability  in  a  Magnetic  Trap  II," 
by  V.  A.  Thomas  and  W,  M,  Nevins,  ERL  Technical  Memorandum  No.  UCB/ERL  M84/46, 
13  June  1984. 

"Vlasov-Poisson  and  Modified  Korteweg-De  Vries  Theory  and  Simulation  of 
Weak  and  Strong  Double  Layers,"  by  K.  Y.  Kim,  ERL  Technical  Memorandum 
No.  UCB/ERL  M84/47,  June  1984. 

"Application  of  Nonlinear  Constants  of  Motion  in  a  Single  Electromagnetic 
Wave  to  the  Study  of  the  Alfv6n-Ion-Cyclotron  Instability,"  by  N.  F.  Otani, 

ERL  Technical  Memorandum  No.  UCB/ERL  M84/49,  16  July  1984. 


Charles  K.  Birdsall 

Professor,  Principal  Investigator 


Thomas  L.  Crystal 
Post-Doctorate 
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SECTION  I:  PLASMA  THEORY  AND  SIMULATION 


A.  Simulation  of  the  Ion-Seam-Driven  Drift  Instability  in  a  Magnetic  Trap 

V.  A  Thomas  (Dr.  W.  M.  Nevins,  LLNL) 


The  following  is  the  abstract  of  a  report  in  preparation. 

Recent  experiments  on  the  TMX-U  device  at  LLNL  have  indicated  the  possi¬ 
bility  of  a  drift  wave  being  driven  unstable  by  the  injection  of  neutral  beams  in 
the  thermal  barrier  region.  A  review  of  linear  theory  is  presented  in  the  local 
approximation.  In  addition,  particle  simulations  are  used  to  understand  the 
nonlinear  characteristics  of  this  instability.  The  particle  simulations  are  per¬ 
formed  using  a  Id  -  3v  electrostatic  particle  simulation  code  and  a  2d  -  3v  elec¬ 
trostatic  particle  simulation  code.  The  instability  is  shown  to  saturate  by  beam 
trapping  nonlinear  electron  effects,  and  nonlinear  motion  in  parallel  with  the 
density  gradient.  Resonant  RrH  motion  is  a  significant  process  for  some  of  the 
beam  particles.  This  process  is  closely  associated  with  the  trapping  of  beam 
ions.  Harmonic  generation  is  also  observed,  as  in  some  of  the  experimented 
data.  Qualitative  nonlinear  theory  is  presented  in  support  of  the  particle  simu¬ 
lations. 


B.  Id  AUvfen-Ioa-Cyclotroa  Instability  Particle  Simulation 

Mels  F.  Otani  (Prof.  C.  K.  Birds  all) 


A  final  report  is  being  readied,  with  title  and  abstract  below. 

Application  of 

Nonlinear  Constants  of  Motion  in  a  Single  electromagnetic  Wave 
to  the  Study  of  the  Alfvfen-lon-CycIotron  Instability 

ABSTRACT 

Constants  of  motion  are  derived  for  particles  moving  in  a  single,  circulary- 
polarized  electromagnetic  wave  of  arbitrary  time-dependence  propagating 
parallel  to  a  uniform  background  magnetic  field.  The  constant  associated  with 
helical  symmetry  is  shown  to  restrict  the  particle  motion  to  a  very  narrow 
region  of  velocity  space.  Features  of  the  slow  time-scale  motion  of  fixed  points 
associated  with  the  existence  of  a  fourth  adiabatic  invariant  are  described  for 
the  case  of  a  slowly  varying  wave.  Characteristics  of  the  particle  motion  thus 
derived  are  applied  to  the  analysis  of  ld-3v  simulations  of  the  saturation  of  the 
Alfvfen-ion-cyclotron  (AIC)  instability  for  a  single  wave.  In  particular,  an  explana¬ 
tion  is  offered  for  the  appearance  of  a  sharp  edge  in  the  velocity  distribution 
function  observed  in  the  simulation. 


C.  Theory  of  Nonmonotonic  Double  Layers 

K  Y  Kim  (Dr.  T.  L.  Crystal) 


The  following  is  an  Abstract  of  a  report  in  preparation. 

We  present  a  simple  graphic  method  of  solving  a  Vlasov-Poisson  system 
associated  with  nonlinear  eigenvalue  conditions  for  arbitrary  potential  struc¬ 
tures.  We  present  a  general  analytic  formulation  for  non-monotonic  double 
layers  and  illustrate  with  some  particular  solutions.  This  class  of  double  layer 
satisfies  the  time  stationary  Vlasov-Poisson  system,  while  requiring  a  Sagdeev 
potential  which  is  a  double-valued  function  of  the  physical  potential:  it  follows 
that  any  distribution  function  having  a  density  representation  as  any  integer  or 
noninteger  power  series  of  potential  can  never  satisfy  the  non-monotonic  double 
layer  boundary  conditions.  A  K-dV  like  equation  is  found  showing  a  relationship 
among  the  speed  of  the  non-monotonic  double  layer,  its  scale  length,  and  its 
degree  of  asymmetry. 


D.  Magnetized  Sheath  Simulations 

Seiji  Ishiguro,  Tetsuo  Katnimura  (IPP,  Nagoya,  Japan) 
(Prof.  C.  K.  BirdsaU) 


Some  rudimentary  simulations  were  done  with  a  magnetized  sheath  at  IPP 
Nagoya.  One  object  is  to  understand  better  the  role  of  the  angle  of  the  magnetic 
field  (with  respect  to  a  bounding  wall)  in  determining  the  sheath  potential  drop 
and  particle  fluxes.  These  runs  were  done  with  our  code  ESI,  modified  to  RAM 
(Reflect-Absorb-Mirror),  done  at  IPP  (Birdsall,  Ohara;  see  our  QPR  3,4  for  1992), 
with  B o  now  added  by  Ishiguro.  Hence,  this  is  a  joint  project,  with  progress 
reported  here. 

The  model  is  shown  in  Figure  1.  Initially,  the  region  0  <  x  <  L/2  is  filled  with 
thermal  electrons  and  ions.  There  is  no  source  so  that  the  initial  plasma  eventu¬ 
ally  decays.  The  wall  at  x  =  L/2  absorbs  all  particles  and  is  floating  (not  con¬ 
nected  to  an  external  circuit);  the  net  surface  charge  ff(L/2)  =  —Es(L/2)  gives 
one  boundary  condition  for  Poisson’s  equation.  The  wall  at  x  =  0  reflects  all  par¬ 
ticles  (changes  sign  of  vs)  and  is  set  at  zero  potential,  <p(0)  =  0,  for  the  second 
boundary  condition. 

RAM  is  obtained  from  ESI  as  follows.  Active  charges  exist  only  for  0  <  x  < 
L/2.  The  charge  density  p  (  0  <  x  <  L/  2)  is  mirrored  across  x  =  L/2  into  L/2  < 
x  <  L  but  not  inverted;  the  potential  is  obtained  for  0  <  x  <  L  as  usual.  The 
active  charges  crossing  x  =  L/2  are  kept  at  x  =  L/2  (and  doubled,  to  account  for 
those  coming  from  the  right),  providing  o(L/2),  which  is  -EX(L/ 2).  The  poten¬ 
tials  at  0  =  x  =  L  are  zero;  volume  charge  weighted  to  x  =  0  is  doubled  to 
account  for  that  weighted  to  x  =  L  (there  is  no  surface  charge  density  at  x  =  0. 
L)- 


The  parameters  used  are: 

L  =  lOOAfl, 

T./Ti  =  100 

rrii/m,  =  100 

v*  =  1 

Upu  —  1 

Uci/ Opi  -  0.  0.1,  0.2,  0.3 
(up/vc)?  =  100,  25,  11.1 
i  =  0.  10° . 80°,  90° 

Some  results  for  an  unmagnetized  ( B0  =  0)  plasma  are  shown  in  Figure  2, 
for  <p(x),  ti*  f  (x),  ^<jPi/t.t,«(x)  intervals  of  u^t  =  100,  and  in  Figure  3  versus 
time.  These  results  are  essentially  just  the  initial  transient.  The  electrons 
charge  the  wall,  dropping  to  roughly  -1/ 2  ln{rr\/ mt)  -  -2.3  in  a  few 

plasma  periods  (i>p,f  »  20).  After  about  time  L/vt ,  =  100,  most  of  the  fastest 
electrons  are  absorbed  and  gradually  rises,  with  oscillations  at  up  and 

(perhaps)  near  u^/2.  The  simulation  would  be  improved  with  a  source,  either 
at  x  =  0  or  distributed  across  part  or  all  of  L.  Nevertheless,  some  results  are  of 
interest.  Note  that  an  edge  in  potential,  density,  and  drift  velocity  propagate 
into  the  plasma  at  sound  speed;  this  is  a  rarefaction  wave,  much  as  seen  in  a 
freely  expanding  plasma.  However,  the  sheath  edge  (where  n*  and  n,  separate) 
is  nearly  stationary,  about  from  the  wall,  as  shown  in  Figure  4. 
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Figure  1.  Reflect-Absorb-Mirror  (RAM)  model,  magnetized.  The  electric 
field  is  Ex(x,t)  only.  The  magnetic  field  is  Bn  independent 

1  L  u 

of  x,  t.  The  wall  at  x  =  ^  (and  at  -  is  absorbing  and 

floating.  This  allows  modeling  the  half  space  as  shown  in 
(b),  rather  than  the  full  space  in  (c),  or  a  full  space 
0  <  x  <  L,  with  a  floating  wall  at  L/2. 


Fl9ure  3-  ^wa 1 1 ’  °wal 1 ’  pnearwall  are  shown  versus  time’  The 

potential  drops  quickly  to  a  value  predicted  by  simple 
sheath  theory,  oscillating  at  about  wpe,  and  modulated 

at  (perhaps)  u.'p./2;  the  slow  rise  reflects  slow  loss  of  the 

fastest  electrons  (there  is  no  source).  The  wall  charges 
negatively  and  the  charge  nearby  shields  the  wall,  as  seen. 


The  magnetized  plasma  ( Bq  *  0)  behaves  similarly,  with  some  notable 
changes.  First,  as  in  Figure  5,  the  ion  drift  velocity  (and  ion  flux)  at  the  absorb¬ 
ing  wall  drops  off,  apparently  as  cos ip,  as  the  magnetization  is  increased. 

Second,  as  in  Figure  6,  (measured  at  a  time  just  past  the  minimum  of 
Vwaii )  appears  rather  insensitive  to  either  ip  or  strength  of  magnetization.  How¬ 
ever,  additional  simulations  showed  that  pvnu  became  positive  with  respect  to 
the  plasma  for  ip  =  90°  ( Bq  parallel  to  the  wall),  as  also  observed  by  L  A. 
Schwager  (Berkeley  PDWl  code,  previous  QPR).  But.  for  ip  =  88°,  ^  is  again 
negative;  the  electrons  simply  take  some  time  to  flow  along  Bo  to  the  wall  The 
time  record  of  ^  now  shows  additional  modulation  apparently  near  uct  >  Up, . 
Hence,  ip  =  90°  appears  singular.  These  90®  effects  (little  or  no  ion  drift  and 
flux,  positive  wall  potential  or  ion  barrier)  certainly  need  more  thorough  investi¬ 
gations.  for  fusion  ion  containment  or  for  device  application  (making  use  of  the 
high  sensitivity  to  ip  for  ip  »  90°). 
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E.  Saturation  Characteristics  of  Counterstreaming  Vann  Electrons 

Any  Wendt  (Prof  C.  K  BirdsaU) 


Counterstreaming  warm  electron  beams  are  unstable  for  ^  lAut,  as  is 
well  known.  Our  study  is  observation  of  growth  as  vt/v0  is  increased  from  zero, 
with  evolution  of  f(t/)  from  two  peaks  toward  (perhaps)  a  Maxwellian.  In  particu¬ 
lar,  f(v  =  0)  is  followed  in  time  and  is  seen  to  overshoot  the  value  for  a  Maxwel¬ 
lian  for  vt/ vq  s  0.3  and  drift  toward  the  value  for  Penrose  stability  for 
vt/v0^  0.5. 

A  report  is  in  preparation. 


F.  Electron  Bernstein  Wave  Investigations  Regarding 
linear  and  Non-linear  Damping 

J.  P.  Lynov  (Ris6)and  Wm.  S  Lawson 


A  modification  of  the  code  PDWl  is  being  used  to  investigate  non-linear 
Bernstein  waves  in  a  plasma.  The  problem  studied  is  that  of  a  spatially  localized 
exciter  emitting  waves  into  a  plasma  close  to  perpendicular  to  the  magnetic 
field.  At  non-linear  amplitudes,  there  will  probably  be  a  wavelength  shift,  and 
reduced  Landau  damping  due  to  the  trapping  of  particles.  It  is  not  presently 
known  at  what  amplitudes  these  effects  will  become  important,  and  it  is  the  pur¬ 
pose  of  this  research  to  find  this  amplitude  as  well  as  the  character  of  these 
effects. 

Also  of  interest,  probably  at  some  future  time,  is  the  mechanism  of  heating 
when  these  waves  do  damp.  This  is  of  possible  importance  to  heating  and 
current  drive  in  Tokamak  plasmas. 


G.  Simulation  of  the  Classical  Pierce  Diode 

T.  L.  Crystal  and  S  Kuhn 


A  report  on  this  project  is  in  preparation,  and  will  be  issued  as  an  ERL 
Memorandum.  This  report  covers  recent  particle  simulations  of  the  "classical" 
Pierce  diode  (1944),  performed  with  the  new  bounded-plasma  code  PDWl.  These 
simulations  (i)  successfully  retrieve  the  linear  behavior  predicted  from  Pierce's 
characteristic  equation,  (ii)  extend  past  the  region  where  linear  instability 
analysis  is  valid  to  recover  the  later  nonlinear  behavior,  and  (iii)  display  the  final 
states  of  the  diode  system  which  might  be  either  time-independent  (steady- 
state)  or  oscillating,  and  which  may  depend  on  the  initial  conditions  given  to  the 
simulation. 
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H.  Simulation  of  the  Plasma-She  ath  and  Presheath  Regions 

Lou  Ann  Schwager  (Prof.  C.  K.  Birds  all) 


Our  particle  simulation  studies  of  the  plasma-sheath  region  compare  favor¬ 
ably  with  analytical  results  of  others.  Our  one-dimensional,  bounded  system  has 
ions  and  electrons  of  equal  temperature  injected  from  one  end  and  an  absorbing 
plate  at  the  other  end,  as  shown  in  Figure  1.  The  code  is  PDWl. 

Equal  fluxes  of  ions  and  electrons  are  injected  at  x=L,  simulating  ions  and 
electrons  escaping  from  a  plasma  center  cell.  Some  electrons  repelled  by  the 
sheath  at  x=0  return  to  the  injection  plane,  where  they  are  thermalized  and  re- 
emitted,  along  with  the  injected  flux.  No  ions  are  returned  to  the  source.  The 
sheath  region  of  a  symmetrical  mirror  reactor  is  modelled;  returning  electrons 
undergo  collisions  in  passing  through  the  center  cell  and  flow  to  the  other  col¬ 
lector  plate. 

At  the  center  cell  plasma  boundary  (injection  plane  in  our  model),  charge 
cannot  accumulate  so  that  the  electric  field  is  made  zero  there.  Thermal  ions 
and  electrons  are  injected  with  equal  densities  and  constant  fluxes  and  tempera¬ 
tures.  The  system  is  initially  empty.  Mass  ratios,  of  40  and  100  are 

used  to  conserve  computation  time. 

The  model  parameters  are: 


number  of  grids.  NG  =  128 
length  L  =  2 

time  step  At  =7.8125  z  10~a 
injected  flux  =  500 
injected  current  =  25 

electron  charge  to  mass  ratio,  (q/m.),  =  -1 
electron  thermal  velocity.  vtt  =  1 

Run  1  Run  2 

rr\/ m,  40  100 

(q/m.)i  0.025  0.010 

vH  0.158  0.100 

number  of  time  steps,  NT  6400  9000 

The  above  parameters  generate  the  following  at  equilibrium. 

Run  1  Run  2 

mi/m,  40  100 

final  number  of  ions  4706  6997 

average  upbt  0.085  0.103 

average  Debye  length,  A/j/Ax  5.90  4.84 

system  length  in  21.7  26.5 


Final  phase  space  and  potential  profiles  are  shown  in  Figures  2  and  3. 

Near  the  emission  plane  the  potential  falls,  accelerating  the  slowest  ions  to 
about  sound  speed,  c,.  as  seen  in  Figures  2b  and  3b.  The  overall  potential  drop 
at  equilibrium  is  1.06  and  1.60  etp/kT  for  the  light  and  heavy  ion  cases,  respec¬ 
tively,  as  observed  in  Figures  2c  and  3c. 

The  above  results  are  compared  with  the  sheath  analysis  of  Emmert,  et  a l 
[l]  for  a  one-dimensional  collisionless  plasma  with  Boltzmann  electrons,  with  a 
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Figure  1.  Bounded  plasma  model. 


PHASE  SPACE  AND  POTENTIAL  PROFILES 
AT  EQUILIBRIUM  FOR  mj/mf=40 


uniform  source.  dependence  on  ffii/m,  for  Z=1  is  obtained  from  their 

equations  and  appears  in  Figure  4  for  Tx/  T,  equal  to  0.10  and  1.0  for  a  source 
width  equal  to  ii  The  of  the  two  PDW1  runs  are  shown  on  the  same  plot, 
even  though  the  source  width  is  2ero.  Also  graphed  is  the  simple  estimate  com¬ 
monly  used,  which  assumes  that  ions  have  a  half-Maxwellian  distribution: 

_  -1.  \M  T* 

KTe  ~  lPn[m  Tx 

(Note:  the  paper  appears  to  have  a  misplaced  exponent.) 

The  paper  also  indicates  that,  in  the  part  of  the  region  where  there  is  no 
plasma  source,  the  potential  profile  is  flat  (except  in  the  sheath).  Hence,  if  one 
chooses  a  narrow  source  width,  say,  10A/j  in  a  system  length  of  100A/>,  the  ana¬ 
lytic  potential  profile  in  Figure  5  is  obtained.  This  result  is  similar  to  the  simula¬ 
tion  profiles  in  Figures  2c  and  3c,  where  the  source  width  is  zero  (planar).  For 
PDW1  with  a  Maxwellian  ion  distribution  emitted  at  the  plasma  boundary,  the 
"point"  source  appears  to  have  a  similar  effective  width.  The  agreement  in  Fig¬ 
ure  4  may  be  fortuitous;  however,  Emmert  et  al.  find  the  overall  potential  drop 
to  be  independent  of  the  source  width  (shape)  (see  beneath  their  Eq.  (35)),  and 
nearly  independent  of  the  system  length  (see  their  Tables  I,  II). 

Future  work  will  add  the  effects  of  secondary  emission  at  the  collection 
plate,  as  well  as  binary  collisions  and  charge  exchange. 


References 

[l]  G.  A.  Emmert,  R.  M.  Wieland,  A.  T.  Mense,  J.  N.  Davidsoa  "Electric  Sheath 
and  Presheath  in  a  Collisionless,  Finite  Ion  Temperature  Plasma."  Phys. 
Fluids  23,  pp.  803-812,  April  1980. 
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Wall  potential  versus  mass  ratio. 
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SECTION  II:  CODE  DEVELOPMENT 


A.  Plasma  Device  Code:  1-d  Vlasov  Simulation  (GASBAG  Code) 

William  S.  Lawson 


A  report  is  nearly  completed,  with  titles  and  abstracts  as  follows. 

A  Method  for  Simulation  of  Vlasov  Systems  with  Discontinuous  Distribution  Functions 

ABSTRACT 

A  method  is  presented  for  simulating  non-periodic  (or  periodic)  Vlasov  sys¬ 
tems  with  discontinuous  distribution  functions.  The  method  is  a  hybrid  of  a 
standard  Vlasov  method  and  the  Waterbag  method;  the  discontinuity  is 
integrated  separately  from  the  distribution  function.  Sample  runs  and  particle 
simulation  comparison  runs  are  included. 

GASBAG  User's  Manual 


Introduction 

This  guide  is  intended  to  be  the  sole  external  documentation  for  the  GAS¬ 
BAG  Vlasov  plasma  simulation  code.  Section  2  is  a  guide  for  the  casual  user  who 
needs  only  enough  information  to  run  the  code.  Section  3  is  for  those  who 
intend  to  use  the  code  as  a  research  tool,  and  will  therefore  want  to  make 
modifications  such  as  new  diagnostics,  and  new  routines  for  loading  or  perturba¬ 
tions.  The  GASBAG  code  was  written  to  be  run  on  the  National  Magnetic  Fusion 
Energy  Compute  r  Center  CRAY  computers.  These  are  vector  machines,  and  the 
code  makes  extensive  use  of  this  vector  capability,  as  well  as  some  software 
unique  to  the  installation.  Anyone  wishing  to  modify  the  code  for  another  instal¬ 
lation  will  find  the  information  in  Section  4  useful. 

It  is  assumed  that  the  reader  understands  the  method  used  by  the  code. 
Using  the  code  without  understanding  the  method  is  likely  to  produce  bad 
results.  For  example,  if  the  string  of  test  particles  finds  its  way  too  far  up  the 
tail  of  the  Vlasov  distribution,  a  small  error  in  the  position  of  the  test  particles 
will  produce  a  large  error  in  the  charge  density,  current,  and  kinetic  energy.  A 
paper  on  the  method  has  been  submitted  for  publication  in  the  Journal  of  Com¬ 
putational  Physics,  and  is  available  from  the  author. 
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B.  PDW1,  Plasma  Device  Code 

William  S.  Lawson 


The  following  are  the  Abstract  and  part  of  the  Introduction  of  the  PDWl 
User's  Manual  now  in  preparation. 


ABSTRACT 

PDWl  is  a  one  dimensional,  magnetized,  electrostatic  particle  simulation 
with  one,  two  or  three  velocity  coordinates,  and  non-periodic  boundary  condi¬ 
tions.  and  includes  a  simple  R,  L.  C,  external  circuit  and  ac/dc  applied  voltage. 
It  is  set  up  to  inject  Maxwellian  distributions  with  arbitrary  drift  velocities  and 
cut-off  velocities,  but  is  easily  modifiable  to  suit  other  distributions.  This  docu¬ 
ment  is  both  a  user's  manual,  and  a  report  on  the  principles  of  operation  of  the 
code.  Inquiries  should  be  directed  to  Prof.  C.  K.  Birdsall. 

Introduction 

This  report  is  intended  as  a  guide  to  the  use  and  modification  of  the  simula¬ 
tion  code  PDWl.  It  includes  detailed  descriptions  of  the  methods  and  algorithms 
used  in  the  code.  Some  familiarity  with  numerical  methods  on  the  part  of  the 
reader  is  assumed.  The  report  is  broken  up  into  six  Sections  and  five  Appen¬ 
dices.  Section  I  is.  of  course,  the  Introduction,  which  describes  the  contents  of 
this  write-up.  and  some  general  programming  concepts  used  in  the  design  of  the 
code.  Section  II  is  a  guide  to  the  mechanics  of  compiling  and  running  the  code. 
Section  III  describes  the  various  parameters  used  to  completely  specify  a  given 
problem.  Section  IV  will  be  helpful  to  anyone  who  wants  to  make  minor 
modifications  to  the  code  (which  will  be  almost  everyone  who  uses  the  code). 
Section  V  is  an  exhaustive  list  of  virtually  all  the  variables  in  the  program.  Sec¬ 
tion  VI  gives  in-depth  descriptions  of  each  of  the  major  subroutines,  including 
the  algorithms,  and  descriptions  of  any  subroutines  called  by  the  routines.  The 
final  five  Sections  are  Appendices.  The  first  Appendix  contains  flow  charts  for  all 
the  major  routines,  the  second  contains  a  description  of  the  workings  of  the  ZED 
interface  (written  by  Niels  Otani),  the  third  contains  some  integrals  which  will 
make  parameter  selection  easier,  the  fourth  contains  documentation  on  a  minor 
but  useful  modification  of  PDWl  named  PDWMAX,  and  the  fifth  contains  some 
examples  which  should  illustrate  some  of  the  problems  PDWl  is  capable  of  solv¬ 
ing. 

PDWl  is  an  acronym  for  Plasma  Device  Workshop.  The  code  was  written  dur¬ 
ing,  and  for  use  in,  a  seminar  workshop  at  the  University  of  California  at  Berke¬ 
ley  on  axially  bounded  plasma  systems  during  the  1983  spring  quarter.  Since 
the  code  is  in  large  measure  a  product  of  the  information  exchanged  during  this 
workshop,  it  was  deemed  appropriate  to  name  it  after  the  workshop. 
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SECTION  ffl:  SUMMARY  OF  JOURNAL  ARTICLES,  REPORTS,  VISITORS.  TALKS 


Journal  Articles 

(1)  K  Y.  Kim,  "Weak  Monotonic  Double  Layers,"  Phys  Letters  Vol.  97A,  No.  1,  2, 
8  August  1983. 

(2)  Yu-Jiuan  Chen,  William  M.  Nevins,  and  Charles  K.  Birdsall,  "Stabilization  of 
the  Lower-Hybrid  Drift  Instability  by  Resonant  Electrons,"  Phys.  Fluids  26. 
September  1983. 


Reports 

(1)  Stephane  Rousset,  "Time-Dependent  Chiid-Langmuir  Diode  Simulation," 
Memo  No.  UCB/ERL  M83/39,  11  July  1983. 

(2)  K.  Y.  Kim,  "Theory  of  Asymmetric  Double  Layers,"  Memo  No.  UCB/ERL 
M83/45,  22  July  1983. 

(3)  S.  Kuhn,  "Linear  Longitudinal  Oscillations  in  Collisionless  Plasma  Diodes 
With  Thin  Sheaths,  Part  I.  Method,"  Memo  No.  UCB/ERL  M83/52,  31  August 
1983. 

(4)  S.  Kuhn,  "General  Path-Integral  Successive-Collision  Solution  of  the 
Bounded  Dynamic  Multi-Swarm  Problem,"  Memo  No.  UCB/ERL  M83/57,  23 
September  1983. 

(5)  S.  Kuhn,  "Linear  Loiigitudinal  Oscillations  in  Collisionless  Plasma  Diodes  with 
Thin  Sheaths,  Part  II.  Application  to  an  Extended  Pierce-Type  Problem." 
Memo  No.  UCB/ERL  M83/81. 


Visitors 

(1)  Dr.  Siegbert  Kuhn  -  from  Professor  F.  Cap's  group  at  the  Institute  for 
Theoretical  Physics,  University  of  Innsbruck,  arrived  in  December  19B2.  He 
is  a  guest  researcher  in  our  group  for  his  sabbatical  year  (1983).  He  is 
working  on  formalisms  for  handling  finite,  bounded  systems  and  aiding  us  in 
our  efforts  at  simulating  axially-bounded  systems  and  sheath  phenomena. 

(2)  Dr.  Jens-Peter  Lynov  -  from  Risp  National  Laboratory,  Roskilde,  Denmark,  to 
work  on  a  variety  of  subjects,  including  heating  and  current  drive  using  off- 
perpendicular  electron  Bernstein  waves  (six  weeks). 


(3)  Dr.  Michael  J.  Gerver  was  here  for  two  weeks  in  July. 


Poster  Papers  for  the  APS  Div.  of  Plasma  Physics  Meeting.  November  1983,  Los 

Angeles 

(1)  (3Z14)  Niels  F.  Otani,  "Simulation  of  the  Alfven  3on  Cyclotron  Instability  in 
Tandem  Mirror  Plasmas." 

(2)  (5W11)  S.  Kuhn.  "The  General  Perturbation  Problem  For  the  One  Dimen¬ 
sional  Collisionless  Plasma  Diode.” 

(3)  (5W12)  C.  K.  Birdsall,  T.  L  Crystal,  S.  Kuhn,  A.  B.  Langdon,  W.  S.  Lawson,  and 
N.  F.  Otani.  "Plasma  Diode  Simulation;  the  PDW1  Code.” 

(4)  (5W13)  W.  S.  Lawson,  "Vlasov  Simulation  on  a  Finite  Grid  With  a  Discontinu¬ 
ous  Distribution  Function." 

(5)  (5W14)  N.  F.  Otani,  "Scaling  of  Magnetized  Double  Layers." 

(6)  (5W15)  K.  Y.  Kim,  "Weak  Asymmetric  Double  Layers  (ADL)." 

(7)  (5W16)  T.  L  Crystal,  S.  Kuhn,  W.  S.  Lawson,  "Pierce  Diode  Instabilities  and 
Extensions;  Simulations  With  PDWl." 

(8)  (7Q6)  Vincent  A.  Thomas.  U.  C.  Berkeley,  William  McCay  Nevins  and  T.  A. 
Casper,  Lawrence  Livermore  National  Laboratory,  ’Ton  Beam  Driven  Drift 
Instability  in  Tandem  Mirrors." 


The  abstracts  are  on  the  page  following. 
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'wort  auppodod  by  DOE  Coo  tract  No.  DE-ATQ3-74ET53064,  comp* 
taaloko  fadUdak  gnmdod  by  NMFECC. 


THE  GENERAL  PERTURBATION  PROBLEM  FOR 
THE  ONE-DIMENSIONAL  COLLISIONLESS  PLASMA  DIODE  • 

S  KUHN.  ERL..  Vrur.  of  Calif..  Berkeley.  CM  94730.  - 
A  method  ij  proposed  for  ireatint  linear  oscillations  in  collision  tern 
one-dimensional  plasma  diodes  with  non-ui.tform  equilibria,  allowing 
for  very  general  inilial.  boundary,  and  eaternal -circuit  conditions.  The 
method,  which  is  an  extension  of  previous  work1-2,  is  baaed  on  the 
basic  (Vlasov  -  Poisson)  and  auxiliary  equations  in  integral  form.  It  is 
potentially  applicable  to  a  variety  of  neutral  and  non-neutral  bounded 
plasma  systems,  such  as  0  machines  or  electron  diodes  The  general 
formalism  is  developed  and  applied  to  the  special  case  of  a  positively 
biased  single  ended  Q  machine. 

*  Work  supported  under  ONR  contract  NOOOI4-77-C-057I,  and  Aus¬ 
trian  Research  Funds  contract  S- 1 8/03 

I  S.  Kuhn.  Proc  19S2  ICPP.  Gdtebotg.  p  1  IP-11-17 
2.  S  Kuhn,  submitted  for  publication. 


PLASMA  DIODE  SIMULATION:  THE  PDW1  CODE 
C  K  BIRDSALL,  T.L.  CRYSTAL.  S.  KUHN.  A.B.  LANGDON,  W.S. 
LAWSON.  N.F  OTANI.  E.R.L..  Uni*  of  Calif..  Berkeley.  CA  94730 .' 
PDWI  solves  for  electron  and  ion  motions  in  the  self-consistent  elec¬ 
trostatic  Acid  of  a  ld(l,2,Jv)  plasma  diode  having  external  series  RLC 
circuit  as  shown  below  Among  many  options,  one  may  model  constant 
current  drive,  one  electrode  Itoating.  or  Axed  potential  drive.  Standard 
code  options  include  various  particle  injection  schemes  (e.g..  drifting 
cutoff1  Maxwellians)  from  either  or  both  electrodes,  and  a  static  mag¬ 
netic  Aeld  B.  at  angle  4  to  the  system  axis.  Recent  applications,  as 
documented  in  related  posters  in  this  session,  include  the  Pierce.  Lang¬ 
muir  (electron)  and  multispecies  diodes;  weak  and  strong  double  layers 
both  with  and  without  B,.  and  magnetic  sheaths  and  Arst  walls.  PDWI 
came  from  a  Plasma  Diode  Workshop,  chief  programmer  W.S.  Lawson. 


*  Work  supported  by  DOE  contract  DE-AT03-76ET53064  and  ONR 
contract  N000I4- 77.00578  Computations  were  done  st  NMFECC, 
Livermore. 


VLASOV  SIMULATION  ON  A  FINITIJ  GRID  WITH 
A  DISCONTINUOUS  DISTRIBUTION  FUNCTION. 

Wat.  S.  Lawson  Electronics  Research  Laboratory.  U.C.  Berkeley.  94720. 
Ordinary  Vlasov  simulation  methods  are  unable  to  deal  with 
non-periodic  boundary  conditions  because  steep  gradients  in  phase 
space  and  discontinuities  are  inherent  in  such  problems  A  method 
which  shows  promise  of  overcoming  most  of  these  difficulties  is  to  ex¬ 
tend  the  distribution  function  post  any  discontinuity,  then  keep  track  of 
the  position  of  the  discontinuity  with  a  string  of  test  particles  Only  the 
part  of  the  distribution  function  which  is  to  one  side  of  the  string  of 
test  particles  is  then  integrated  in  calculating  the  charge  density  The 
unphysical  extension  of  the  distribution  function  does  not  influence  the 
physics,  yet  it  allows  Anne  difference  methods  to  be  used  to  integrate 
the  Vlasov  equation  Well-behaved  Vlasov  integrators  have  been 
found,  and  results  from  a  code  based  on  this  method  (GASBAG)  com- 
pare  very  well  with  particle  simulation  results  While  the  accuracy  of 
the  Vlasov  simulation  is  not  better  than  that  of  the  particles,  the  noise 
level  is  far  lower 

•  Work  performed  under  Dept  of  Energy  contract  DE-AT03- 
74ET53064  using  the  NMFECC  computers 


SCALING  OF  MAGNETIZED  DOUBLE  LAYERS. 

N.  OTANI.  E.B.L..  Unix  of  California  Berkeley.  CA  94720  • 

Preliminary  reaults  from  one-dimensional  bounded  electrostatic  panicle 
aunuiauoas  employing  constant  current  boundary  conditions  ere  dis¬ 
cussed.  We  obeerve  magneuud  double  layers  persisting  for  several 
hundred  electron  plasma  periods  A  large  current-driven  positive  poten 
Us!  structure  eventually  develops  near  the  emitting  wall,  limiting  the 
useful  duration  of  our  simulations  Obliquely  magnetised  double  layers 
scale  spatially  as  the  Debye  ler«lh  during  the  Arst  few  ion -cyclotron 
periods. 

•  Work  supported  by  DOE  contract  DE-AT03-76ET53064, 
computational  facilities  provided  by  NMFECC,  Livermore 


WEAR  ASYMMETRIC  DOUBLE  LAYERS  (ADL) 
K.  Y.  KIM,  £J(X.  Uatr.  tfCakf.,  Berkeley  CA  94730.' 

A.  Otar  analytic  formulation  for  ADL  involves  a  third  constant  of 
erosion.  agnCr-a.)  taken  x.  locates  Iks  ADL  minimum  (or  max¬ 
imum)  potential,  in  addition  to  tbs  usual  constants,  total -energy  and 
sga(r).  (My  ik  this  way  can  wa  guarantee  ikat  (I)  Ike  distribution 
AiUlgS  — '"**■  tbs  Vlasov  eqo  and  ie  continuous  throughout  phase 
mace,  and  (2)  the  Sgitffgv  potential  ie  double  valued.  Using  amplitude 
dependent  (nonlinear)  boundary  wwAikw*  with  these  distribution 
tactions  yields  the  determinants!  aquations  producing  the  ADL. 

R  pardda  simulations  of  ADL  have  been  successful  in  short  systems 
(i.  <  gOJtp)  and  with  low  drift  vtioauca  <  V,  <  0  5  P,  for  the  elec¬ 
trons)  Wa  have  simulation  reaults  for  systems  driven  by  constant 
current  and  by  constant  applied  voltage,  both  with  sod  without  applied 
magnetic  Bold. 

*  Work  potformod  under  Dept  of  Energy  contract  DE-AT03- 
7MT5JOM  using  NMFBCC  computed  at  Livermore 


PIERCE  DIODE  INSTABILITIES  AND  EXTENSIONS: 
SIMULATIONS  WITH  PDWI  •  T.L  CRYSTAL.  S  KUHN*.  Wm 
LAWSON.  ERL.  V.  C  Berkeley.  CA  94720 

Particle  annuls ootu  of  the  clawcsl  (xhwt-urtuiicdl  Pierce  diode  recov¬ 
er  linear  eigenmodae  and  agenfrequcnciei  as  predicted  by  linear  theory 
Stability  of  this  special  plasma  diode  is  governed  by  die  free  parameter 
a  “w, reflecting  effects  of  both  the  ylasma  properties  and  the 
boundary  conditions  When  the  diode  is  unstable,  the  system  evolves 
to  ■  nonlinear  Ana!  state  which  can  depend  fundamentally  on  the  motel 
conditions  Simulations  including  external  RLC  circuit  elements 
demonstrate  that  these  seriously  modify  both  the  character  and  the  sta¬ 
bility  of  the  diode  system's  response. 


•  Work  supported  under  DOE  contract  DE-AT03-76ET53064.  and 
ONR  contract  N 000 U-77-C -0578  One  of  the  authors  (SKI  ack¬ 
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